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U ' Abstract 

■B 

Two related open problems in the theory of iD Navier-Stokes turbulence 
are discussed in this paper. The first is the phenomenon of intermittency 
in the dissipation field. Dissipation-range intermittency was first discovered 
^ ' experimentally by Batchclor and Townsend over fifty years ago. ft is charac- 

\^ , terized by spatio-temporal binary behaviour in which long, quiescent periods 

\l ' in the velocity signal are interrupted by short, active 'events' during which 

. _^ , there are violent fluctuations away from the average. The second and related 

^— V ' problem is whether solutions of the 3D Navier-Stokes equations develop fi- 

^+ , nite time singularities during these events. This paper shows that Leray's 

^^ ' weak solutions of the three-dimensional incompressible Navier-Stokes equa- 

^^ , tions can have a binary character in time. The time-axis is split into 'good' 

and 'bad' intervals: on the 'good' intervals solutions are bounded and reg- 
ular, whereas singularities are still possible within the 'bad' intervals. An 
estimate for the width of the latter is very small and decreases with increas- 
ing Reynolds number, ft also decreases relative to the lengths of the good 
intervals as the Reynolds number increases. Within these 'bad' intervals, 
lower bounds on the local energy dissipation rate and other quantities, such 
C^ ' as ||m(-, t)\\oo and ||Vm(-, i)||oo, are very large, resulting in strong dynam- 

ics at sub-Kolmogorov scales. Intersections of bad intervals for n > 1 are 
related to Scheffer's potentially singular set in time. It is also proved that 
the Navier-Stokes equations are conditionally regular provided, in a given 
'bad' interval; the energy has a lower bound that is decaying exponentially 
in time. 
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Fig. 1. A typical example of dissipation-range intermittency from wind tunnel 
turbulence where hot wire anemometry has been used to measure the longitudinal 
velocity derivative at a single point (D. Hurst and J. C. Vassilicos). The horizontal 
axis spans 8 integral time scales. The Taylor micro-scale based Reynolds number 
is about 200. 



1. Introduction 



The questions to be addressed in this paper concern the nature and 
behaviour of intermittent high Reynolds number solutions of the three- 
dimensional Navier-Stokes equations. Relatively quiescent flows can exist 
in nature at high Reynolds numbers. There are also flows that appear tur- 
bulent for all practical purposes but are nevertheless smooth at appropri- 
ately small length scales. Typically turbulent high Reynolds number Navier- 
Stokes flows, however, generally display a specific hallmark which is called 
dissipation-range intermittency. This was first discovered by Batchelor and 
Townsend (Jj and manifests itself in violent fluctuations of very short dura- 
tion in the energy dissipation rate. These fluctuations away from the aver- 
age are interspersed by quieter, longer periods in the dynamics. The data 
in Figure Q is an illustration of a typically intermittent signal represent- 
ing a velocity derivative versus time recorded at a single point in space. 
For theoretical studies of the Navier-Stokes equations it is usual to express 
energy dissipation in the L^-volume-integrated sense. Postponing the full 
definition of system variables, if the global energy Hq and enstrophy Hi for 
a three-dimensional Navier-Stokes velocity field u{x, t) are defined as 



Hoit) 



\uVdV 



Hi[t) = / \Vu\^dV 
Jv 



(1) 



then it is well known that wliile the energy Hq is a uniformly bounded 
function of time, only the long-time averaged energy dissipation rate eav = 
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i/L'^ (Hi) is known to be bounded |2], whereas the behaviour of e{t) ~ 
vL~^Hi{t) pointwise in time may be wildly fluctuating or even singular. 
The long-time average (■) is defined later in this section. 

Whether Hi{t) becomes singular in a finite time is an open question 
intimately related to the Navier-Stokes regularity problem. So long as Hi{t) 
is finite the solution is smooth and unique; any finite time singularity must 
be accompanied by a divergence of Hi[t) at that time. From Leray onwards 
|2], this question has led to a long and rich literature on the nature of 
weak solutions [5| | ll l il lM7ll51l^l(Jllll| . In terms of physical length scales, the 
boundedncss of {Hi) allows the time averaged turbulent energy dissipation 
rate per unit volume €av = vL~^ {Hi) to be used naturally in forming the 
inverse Kolmogorov length 

Spiky behaviour in Hi (i) , causing loss of resolution in large-scale computa- 
tions, could mean that significant energy lies in wave-numbers k > rj^ in 
the dissipation range of the energy spectrum. The ubiquity of dissipation- 
range intermittency in turbulent flows suggests that it should occur natu- 
rally in mathematical analyses of the Navier-Stokes equations. Strong tem- 
poral excursions in Hi{t) are clear candidates for the formation of singu- 
larities and may be related to potentially singular solutions of the three- 
dimensional Euler equations |12II13II14II11) . although nothing has been rig- 
orously proved in this respect. Beginning with Leray 's seminal paper in 
1934, rigorous methods of analysis on the full three-dimensional domain 
have led to seventy years of literature on the Navier-Stokes regularity prob- 
lem but have yet to settle this question definitively. Short-time regularity 
has also been known for many years, as have various interesting partial and 

Batchelor and Townsend J] also suggested that energy dissipation is 
not distributed evenly across the full three-dimensional spatial domain but 
is clustered into smaller spots in the flow with the energy associated with 
the small-scale components being distributed unevenly in space and roughly 
confined to regions which become smaller with eddy size |18| . In contrast 
to Kolmogorov's traditional theory that implies that energy dissipation is 
space-filling ^Hl, Mandelbrot suggested that the spatial set on which it 
occurs is actually fractal [501 . In experimental investigations of the energy 
dissipation rate in several laboratory flows, and in the atmospheric surface 
layer, Sreenivasan & Meneveau [211122) interpreted the evident intermittent 
nature of their signals in terms of multifractals (see also [221 )• Zeff et al [2^ 
have shown how more recent technical advances have made it possible to 
measure each derivative of all three velocity components to obtain a fuller 
experimental picture of the energy dissipation at a point in a flow. In general, 
numerical simulations and experiments suggest that respectively quasi-one 
and two-dimensional tubes and sheets are the favoured low-dimensional sets 
on which vorticity and strain appear to accumulate [5H1I26II27) although, 
as Galanti & Tsinober |2H1 and Tsinober [201 have pointed out, there are 
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significant differences between ttiese two sets. It is also true tlrat explaining 
these phenomena in the simple geometrical terms of tubes and sheets is 
a visual over-simplification of much more complicated dynamical spatial 
structures at small scales. Examples of this are the spiral vortex structures 
introduced by Lundgren j^Ol and discussed in detail by Vassilicos & Hunt 
|31j . Flohr & Vassilicos [22 and Angilella & Vassilicos 



Dissipation-range intermittency is a well established, experimentally ob- 
servable phenomenon; its appearance in systems other than the Navier- 
Stokes equations has been discussed in an early and easily accessible paper 
by Frisch & Morf [2H1- One symptom of its occurrence is the deviation of 
the 'flatness' of a velocity signal (the ratio of the 4th order moment to 
the square of the second order moment) from the value of 3 that usually 
holds for Gaussian statistics |21] . More subtle is the phenomenon of inertial- 
rangc intermittency that has exercised the ingenuity of that section of the 
physics community that focuses on scaling methods. In the inertial range, 
Kolmogorov's theory predicts that the exponent, ^p, of the pth velocity 
structure function should vary linearly with p, whereas experimental data 
shows that the ((^p, p) relation is a concave curve lying below the line p/3 for 
p > 3. This departure from Kolmogorov scaling in the inertial range, and 
therefore from the five-thirds law, is termed inertial-range intermittency. An 
extensive literature is quoted in Frisch's book |19| . Studies in weak turbu- 
lence, applicable to predominantly dispersive systems, have been pioneered 
by Zakharov. These ideas can be found in papers by Zakharov, L'vov and 
Falkovich PTj and Zakharov |38j . 

To prove that solutions of the Navier-Stokes equations (usually taken 
in a periodic box |51l^lOII39p are typically intermittent in space-time poses 
formidable technical challenges to the mathematician. Analysis on time- 
evolving fractal domains with ill-defined boundary conditions is not ad- 
vanced enough to gain rigorous results by concentrating on one fractal 'spot' 
in the flow. The partial regularity result of SchefFer ^^l proving that the 
potentially singular set in time alone has zero half-dimensional HausdorfF 
measure has been of great influence (see also JBl for a short proof) . Follow- 
ing this, the most significant space-time result has been that of Caffarelli, 
Kohn & Nirenberg ^21 who showed that the potentially singular set in 
space-time has zero one-dimensional Hausdorff measure. This implies that 
if singularities do exist they must be relatively rare. Lin |40| and Choe & 
Lewis |41| have recently provided shorter proofs of this result. 

The more realistic option adopted here is to show that solutions of the 
Navier-Stokes equations can have a binary nature in time in which the time- 
axis is divided into what are designated as good and bad intervals. On the 
good intervals the Navier-Stokes equations are uniformly regular. The bad 
intervals are shown to be very small in width with an upper bound that 
decreases with increasing Reynolds number and which also decreases rela- 
tive to the widths of the good intervals. Within the bad intervals very large 
lower bounds are shown to exist on both the local-in-time energy dissipation 
rate and several other quantities, such as ||m(-, i)||oo and [[Vtt(-, t)||oo- The 
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corresponding local length scales within these intervals are, at best, compar- 
atively much smaller than the Kolmogorov length. The regularity question 
within the bad intervals is still open, so only weak solutions are known 
to exist there. The great difficulties encountered by computational fluid 
dynamicists in resolving turbulent flows even for modestly high Reynolds 
numbers could be because of this binary behaviour. 

These results, which are summarized in §2, have been obtained through 
the use of a set of quantities K„(i) that have been introduced in previous pa- 
pers |59,42,43 . The more physical aspects of these ideas have been laid out 
previously in a short paper |44j ; the present paper gives a detailed and more 
advanced account of the methods and results reported there. The K„(t) have 
the dimensions of inverse lengths and are formed from ratios of L^-norms of 
derivatives of the velocity field. Together with the periodicity of the domain, 
the X^-spatial integration within the k„ (i) means that spatially intermittent 
effects are included implicitly and cannot be averaged away. Clearly they 
are not the same quantities as those measured by experimentalists, such as 
the energy dissipation rate, but estimates for them are rigorous, making no 
appeal to any approximations, and ultimately lead to information on the 
energy dissipation. Physically they can be considered as a measure of the 
2rith moment of the energy spectrum. Their time-dependence is explicit and 
their long-time averages (k„) are uniformly bounded |43| . Pointwise in time 
their binary nature appears for each value of n > 2. 

The phrase 'solutions can have a binary nature in time' has been used 
above in the following sense: if no bad intervals occur, then solutions of the 
Navier-Stokes equations are bounded for all time. In this sense, the results 
in this paper are different from conventional short-time regularity proofs be- 
cause loss of regularity can only occur in the bad intervals. This is consistent 
with Scheffer's partial regularity result |15j : if singularities exist at points 
in time then these must be clustered within the intersection of the bad 
intervals. They are also consistent with the well known problems of compu- 
tational resolution in three-dimensional turbulence: the bounds indicate a 
structure so fine that it would be extremely difficult to distinguish between 
regular and singular solutions, despite the possibility of a lack of sharp- 
ness. At this point it is appropriate to describe the Navier-Stokes system of 



Table 1. Definitions of the main parameters in the paper. 
Quantity Definition Comment 



Box length 


L 


Forcing length scale 


I 


Average forcing 


p-L-^mi 


Narrow-band forcing 


\\T\\2 ~ 1^ II V JII2 


Average velocity 


U'=L-m\u\\l) 


Grashof No 


Gr = /^3i/-2 


Reynolds No gS] 


Re = Ulv-^ 



< L/2-K 



Gr^/^ <cReas Gr 
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partial differential equations considered on a periodic cube V — [0, L]'^ 

Ut + u-\7u = iyAu-\7p+ f{x), V-m = (3) 

where v is the kinematic viscosity and p is the pressure. The applied body 
force f{x) is taken to be mean-zero and divergence-free so, without loss 
of generality, the solution u{x,t) is mean-zero at all times. For simplicity 
narrow-band body forces with a single length scale i are considered; that is, 
with Fourier components only at wave-number k = £~^ and i < L/2i:. For 
finite energy initial data the Navier-Stokes equations admit weak solutions 
in L'^{V) at each instant of time, with finite time integrals of the L^-norms of 
the velocity gradients. With the assumption of narrow-band forcing, norms 
of gradients oi f{x) are related to the norm oi f{x); these can all be found in 
Tabled Also found there are the definitions of the root-mean-square velocity 
scale U and the Reynolds number Re. The angled brackets (■) denote the 
long-time average 

(<?>(■))= LIM,^^ ^ <?>(s)ds). (4) 

LIM is a generalized long-time limit for functionals of (weak) statistical solu- 
tions of the Navier-Stokes equations i5i|6j. The square of the L^-norm || ■ HI 
is defined as 

Il/ll2= / \f\'dV (5) 

Jv 

The Grashof number Gr is the natural control parameter, not the Reynolds 

number Re, but it is clear that high Reynolds number solutions may be 

achieved if the Grashof number Gr is sufficiently high; indeed, Doering & 

Foias 03 have proved that for body-forced Navier-Stokes flows such as these 

Gr ' < cRe as Gr -^ oo. Next are defined the time dependent quantities 

Hn{t)-- I |V"M(a;,t)pdF = Vfc^"|M(fc,t)p. (6) 

where u{k,t) is the Fourier transform of u{x,t). For higher derivatives of 
u and /, the important quantities that will be used in this paper are 

i^n== / (|V"Mp + r2|V"/P) rf^ (7) 



where r is a characteristic time the origin of which is discussed later in ^ 
In approaching this problem conventional L^-norms have been used in (O 
to avoid difficulties with the pressure field. 

Not all interesting small-scale spatial behaviour can be averaged away 
on a finite periodic box; spatial events must show up in some temporal 
manner. The next step is to define the quantities 

^nAt) = (i^Ji^,.)'/'^"-'^^ (8) 
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The particular one of interest is k„ ~ k„^o where the r ~ label has been 
dropped for convenience in the first, k^" can be interpreted as being related 
to the 2nth moment of the energy spectrum )48| : this can be seen from 
writing jvp ~ \u\'^ + |/p and then applying Parseval's equality to (jHJ (with 

Moreover, the k„ are ordered in magnitude for all t > 

L^^ < Kl ■ ■ ■ < Kn < Kn+l < ... (10) 

which is simply a result of Holder's inequality. The full k„ ,, are also ordered 
such that Kn.r < n-n.r+i for r + 1 < n. 



Table 2. Definitions of the main variables and constants in the paper. The pa- 
rameter S, which lies in the range < (5 < i, is ignored hereafter. 

Definition 
Hn H„{t) = Jy\V"u\UV n>0 

Eai) tav = I'L (Hi) 

F„ F„(i) = L (|V"u|2 + r^lWI') dV 

uoT = Gr-(i/^+^) < 5 < i 

^u Kn(t) = (i^n/i^o)''"" 

A„ -^1 = 3 - 2;7 + - 



2. Summary of the main results of the paper 

Sections 3-5 of this paper give a full account of the ideas with proofs from 
first principles. This section has been introduced for readers who prefer 
to peruse the proofs later. The next sub-section on long-time averages is 
followed by a second describing where problems with regularity lie. The 
third describes how the time-axis can be split into two types of interval, 
'good' and 'bad'. ii2. 41 summarises the dynamics on the bad intervals and 
the very large lower bounds that can be found within these. 

2.1. Long-time averages 

The best known long-time average in Navier-Stokes analysis is the ex- 
plicit upper bound on the time averaged dissipation rate eav = vL~^ {Hi) 
found from Leray's energy inequality [21 

hHo<-yHi+Hl'^\\f\\2 (11) 



8 J. D. Gibbon and Charles R. Doering 

Using Doering and Foias's result 0^1 that Gr < cRe '-, it is easily shown 
that as Gr — > oo 

(Hi) <ciy^L^r'^Rc^ (12) 

Recall that the long-time average (•) is defined in Q). The above Re"^ long- 
time average estimate for (Hi) leads to an upper bound on the energy 
dissipation rate eav < c u^l^^Ke and thenee to the conventional estimate 
on the inverse Kolmogorov length rj^^ defined in jSJ 

77^i=(^)'^'<riRe3/4 (13) 

Foias, Guillope & Temam's generalization of (|12|l . when the inequality Gr < 
cRe^'^ is apphed |7j, can be expressed as (see Theorem |21 in H3.3|) 

i /f^\ < Cniir'^fi^^^Re^ (14) 

When n = 1, p4|) recovers the sharp result of Doering and Foias jl^ for 
(Fi), except for a spurious volume factor (Li^^)^ on the right hand side. 

A closely connected and important result used in this paper is the bound- 
edness of the long-time averages of the k„, for n > 1 taken with i = L/2Tr 
(see gS] and Theorem [T] in ^J^ 

(LHin) < c„Re^" A„ = 3 - ;^ + - (15) 

Zn n 

where (5 is a small parameter lying in the range < S < i. From either 
(fT^ or (fT5|l an upper bound can also found on (||w||oo); see Theorem El in 
i)3.3l and Tabled Note that the case n = 1 gives (Lui) < ciRc ' , which is 
consistent with the traditional scaling of Re ' for the Taylor micro-scale. 
Boundedness of the long-time average of k„ does not, of course, imply that 
the Kn are bounded pointwise in time. 

2.2. Problems with regularity: an illustration 

Normal practice has been to consider the time evolution of the F^ us- 
ing differential inequalities. With variations, this has been standard the 
approach taken since the early days of the subject [5IIUlllOII39| . One such 
inequality is 

iF„ < -ii^F„+i + (c„i/-1||m||^ + i^r^Rc) F„ (16) 

where ||m(-, t)\\oo is, in effect, the peak velocity on the whole domain. The 
reader can find the precise derivation of (|16|l in Proposition ^ in i)3.2l The 
right hand side has two dominant terms; one negative associated with the 
dissipation, and the dominant positive ||it||;;Q-term. Rewriting the -Fn+i term 
in 1)16(1 in terms of the k„ defined in (jHJ 

F^+.^nir^] F^. (17) 
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and using a Sobolev inequality ||m||^ < ck^Fq for n>2, (|16|l becomes 

r / \2(n+l) ^ 

hFn < < -hy>il ( -^ j + ^nV-^^lF^ + Vt-''^^ \ Fn- (18) 

From pU|) the ratio k„+i/k„ has a lower bound of unity, k„+i/k„ > 1, for 
all periodic divergence-free functions. This reduces (|18|l to 



Fn < {-h^^^l + CnV'^KlF^ + vr'-Ke) F„. (19) 



The manifestly negative term ~ k^ in H19|l is not sufficient to control the 
K^ term: arbitrarily large initial data on F„ can be chosen that makes the 
right hand side positive. Despite the finiteness of the time averages (|14|) and 
(tT5|l . this leads to a failure to control either F„ or k„, other than for short 
times or small initial data [3]. 

2.3. The potentially binary nature of the time-axis 

One way of proceeding with this difficult problem is to ask whether 
the lower bound of unity on k„_|_i/«;„ could be raised, thereby effectively 
increasing the dissipation in 118|l . Batchelor & Townsend PP experimentally 
identified a similar quantity to k„_|_i/k„ for n = 2, 3 that was larger than 
expected for Gaussian data. The principal result of this paper, proved in 21 
is that an improved lower bound on k„_|_i/k„ can be found which is valid 
only on sections of the time-axis. The estimates in Theorcm|21in 21 the key 
result of the paper, show that for Navier-Stokes weak solutions 



Kn+l 



{LKnTRc-^" ) > (20) 



where the real parameter fi can take any value in the range < /i < 1. 
The Cn are the same as in (|15|) . Given that the long-time average in poil 
is non-negative means that there must be intervals of the time-axis, called 
good intervals, where the inequality 

Kn+l \ ^ fr ., \Mt3„-A 



>(LK„)^Re-^" (21) 

holds. It is easily seen that when (|21|l is applied to H18|) on these intervals 
the strength of the dissipation is increased. This applies at small scales 
{Lkji > c„Re "''') where the lower bound on k„+i/k„ in H21|) is raised away 
from unity. The divisor (Fq) within k„ is bounded both above and below. 
Thus H18|l can be turned into a proper differential inequality in F„; the 
reader can turn to t|4.1l for details. The result, which is intuitively obvious 
from (|18|1 . is that the negative dissipation term is stronger than the positive 
nonlinear term when 

^ > ^^T^ (22) 

2(?i+ 1) 
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so no singularities are possible on these intervals (see t l4.1|l . 

The integrand in H2()(l cannot be guaranteed to be always positive, how- 
ever; so intervals on the rest of the time-axis, where it could be negative, 
are designated as bad intervals. On these the reverse inequality must be 
satisfied 



Hn+l 



Kn 



< {Lfinf Re" 



Because its is always true that k„+i > k„, then 

LK„(i) > c„Re^"/^ 



(23) 



(24) 



on these intervals. It is, of course, possible that there are no bad intervals 
and that the whole time-axis is 'good'. The positivity of the average in (|2()|l . 
however, ensures that the complete time-axis cannot be 'bad'. This paper is 
based on the worst-case supposition that bad intervals exist and need to be 
dealt with accordingly. There appears to be no further information within 
(|2()|l as to their distribution, which may depend on the value of n. 



2.4- Dynamics on the bad intervals 



Knit) 



Re 



>.n/l^ 



fle^ 





- 


■ Ik ^ Rp^^/^ / 






\ i 






/ 
/ 


\ 1 

\ i 

Long-time avaragc 




\ 


-^ 




\J 




»- 



Figure 2: A descriptive picture, not to scale, of good/bad intervals for some value 
of n > 2; constants have been omitted. 



Various questions remain. Firstly, if the bad intervals exist, are they finite 
in width? The answer to this is in the affirmative. Estimates for the widths 
of bad intervals {At)b are displayed in Table |5l (see Theorem |5l in t |5.2|l 
with upper bounds on fj, lying in two ranges together with a lower bound 
given in 1)21(1 . It can be seen that these widths are exceedingly small; in fact 
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(Z\i)(, — > as Re ^ oo. The large lower bound on k„ in H24(l indicates the 
predominance of high wave- numbers within these very short intervals. 

Figure [Ql is a descriptive picture of a typical distribution of good/bad 
intervals on the t-axis. It is drawn in such a way that K„(i) looks a relatively 
flat function. While this is probable, some artistic licence has been taken for 
the following reasons. The lower bound in l|24|l makes values of K„(t) much 
larger than the upper bound on the time average p5|l . so it has to spend 
relatively long intervals of time in the good to recompense. Nevertheless, 
this does not prove that it is quiescent in the central part of the interval. 
As is shown in 21 there is enough freedom within the upper bound on k„, 
under the constraint of the long-time average, to allow it to reach large 
enough values so that it can connect with the next bad interval. While we 
believe that it is likely that k„ is flat in the central region, we have been 
unable to prove this. It is pathologically possible that enough fine structure 
might exist within the good interval without violating the average. 

The positions of the intervals may differ as n varies so a new figure is 
needed for each n; Lemmas shows that if any one k„ is bounded then all 
are bounded. Only if the intersection of all bad intervals is itself 'bad' is 
a singularity possible; i)5.^SI addresses this question. The set containing the 
intersection of all the bad intervals for all ri > 1 is designated there as iS*^""-* 
and is related to Scheffer's potentially singular set ^^l- Taking this into 
account narrows the range of /x to | < /x < i . 

The bad intervals (and their intersection) can be divided into sub- intervals; 
for the ith bad interval let us take the jth sub- interval on which Fn > 0. 
This we call the ijth 'dangerous' sub-interval of width (zAt^'')- Singularities 
can occur on any interval within these because Fn is increasing, whereas 
sub-intervals on which Fn < no singularities can occur because Fn is 
bounded (n > 1) by its initial value at the start of the sub-interval. Esti- 
mates for the width of dangerous sub-intervals, of width (Z\i!|l-'), are smaller 
than those for (Z\t)h and can be found in Theorem [7| of i)5.4l fsec Table |3}. 
It is then possible to find very large lower bounds within (At^^^ ) on vari- 
ous quantities such as an equivalent of the inverse Kolmogorov length 77^ , 
based upon lower bounds on Fi within (Z\t!|l-'); it is here where the inter- 
mittency of the dissipation field shows up. Additionally, the peak velocity 
and the peak velocity gradient matrix also have very large lower bounds. 
These are displayed in Tabled 

Conclusions that can be drawn from Tables |3| 0] and ]E\ at the back of 
this section are: 



1. The bad intervals are exceptionally narrow for large Re. 

2. The action within these intervals is intense; the lower bound on Fi within 
them illustrates the strength of the dissipation field there. 

3. The dynamics within these intervals is so fine, even if no singularities 
occur, that they would be exceptionally difficult to resolve numerically. 
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Thirdly Theorem 01 bounds below, in the average sense, the ratio of the 
widths of the good and bad intervals 



lfl^>Re^"(^-i) (25) 

showing that relative to the bad intervals, the widths of the good increase 
with Re. 

Fourthly, the final result of ^Hl displayed in Theorem |H1 is a conditional 
regularity result. Assume that the energy Ho(t) has a lower bound within 
the dangerous sub-intervals, At\.' = t — t^-' , of the form 

Ha{t) > ifo(4'-'')e'"°^'=^*+"' . (26) 

Then solutions of the Navier-Stokes equations are regular there, t^-' is the 
initial time for the dangerous sub-interval. The very large initial conditions 
on the energy ifo(^o"') ^^ ^^^ junction of the intervals is the main obstacle 
to completing a regularity proof. 



Table 3. A comparison between the upper bounds on the long-time averages in 
column two (see Theorems Q and |21 in tl3.8|l and the very large lower bounds in 
dangerous sub-intervals (At^^-') in column three (see Theorem in >I5.4^ . Notice 
there the very large lower bounds on _Fi illustrating the strong intermittency in 
the dissipation field. All multiplicative constants have been omitted to save space. 
In addition £ = L/2n. 

Long-time average (■) On {At^_^-') intervals 

Energy moments (LKn) < Re^" Lk„ > Re*+''" 

Enstrophy L^^ (Fi) < uj'^Re^ L-^Fi > t^gRe'*+''" 

(Kolmogorov scale) "^ Lr?^,^ < Re^/"* Lr]-^ > Re(''+''")/^ 

(Taylor micro-scale)"^ (Lki) < Re^^^ Lm > Re*""/^ 

Peak velocity {\\u\\oo) < LtooRe^ \\u\\^ > LtuoRe'^'^''" 

Vel. gradient matrix /[|Vu||i^^\ < t^oRe^ ||Vu||oo > c^oRe""''''" 



Table 4. Definitions of A„, a„ and b„. On the intersection set 5'°°' of tl5.3l A„ 
can be replaced by A„ . 



Definition Range 



Xn 


A„ = 


a-n. 


On = 


bn 


&„ = 



_ 5 
A. 



"^ 2n, 



+ 1 I 2n-2 \ 
I y2n-lj 



a„> 

bn > an 
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Table 5. Estimates of widths of bad intervals, their sub-intervals {At',^-'), and 
corresponding ranges of a„, b„ &i fi (see Theorem |K| in tl5.2l and Theorem |7| in 
H5.4|l . It is always true that b„ > a„ (see Table |1J. For the case n = 1, ^ > i. 

Range of /x an, b„ Width of interval n > 2 

M<A„+i(^) a„ > 1 LJo{At)t<Re-''" 

A„+i (^^) < ^ < A„+i (^^^) 0<a„<l LJo{At)t<Re'^[nRe 

fo„ > 1 Luo(AtY) < Re-''" 

<M (Zty >(Zty,Re""(M-0 



2(,i+i) ^ ^^ y^'-'g - y^'-'b 



3. Standard estimates 

3.1. The forcing does not dominate the fluid 
The technical parts of this paper revolve around the quantities 

Fn{t)^Hn + r^\\S/-f\\l (27) 

where the /J„ are defined in ^. The F„ contain the fluid velocity derivatives 
and those of the forcing, although the latter has been assumed to have a 
narrow-band character, as shown in Tabled They arc included within Fn 
in order to circumvent problems that may arise when dividing by these 
(squared) semi-norms. Once these terms have been introduced it is necessary 
to demonstrate that they do not dominate the fluid |121 7 a- point that has also 
made by Tsinober |46| in numerical computations. The characteristic time 
T will be chosen for convenience but so long as t 7^ 0, the Fn are bounded 
away from zero by the explicit value t"^ L^ i~'^^ f'^ . Moreover, r may be chosen 
to depend on the parameters of the problem such that (F„) ~ (Hn) as 
Gr ^ 00. To see how to achieve this, let us define 

r = €2^-iGr-(*+i/2) (28) 

with (5 > 0, which is a parameter yet to be determined. Then the additional 
term in H27|l is 



r^ 



IWII2 = i'i^-'€"-'"/'Gr-(2*+i) - i.2^-(2"+2)^3g^i-2i^ (29) 

Now it has also been shown by Doering & Foias 03 that the energy dissi- 
pation rate eav = vL^^ (^1) has a lower bound for high Gr 

eav>cu^r^L-^Gv (30) 

which can be used in the far right hand side of H29|) 

2||v7n j;||2 ^ „ , «-(2n-l) r 4, -lr^„-25 



T2||V"/ll2<cea.r(2"-i)L4i.-iG 



r 



.r-^' 



< c {M-^f"-'\Hn)Gr-^' (31) 
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where Poincarc's inequality has been used at the last step. Hence, for any 
i5 > the additional forcing term in H27() is seen to be negligible with 
respect to (i/„) as Gr -^ oo. The parameter 5 is left arbitrary at this stage, 
although it will be restricted further in the course of proving results in the 
next section. Our interest lies in results for high Gr so correction terms 
described above will be ignored and it can safely be said that for 6 > 

{Fi) <ciy^L^e-^Re^ (32) 

This is Leray's result for weak solutions with narrow-band forcing included 
in a rational manner; the next section deals with long-time averages of other 
quantities. 

3.2. The Fn-ladder 

In the calculations that follow the Hn are formally manipulated even 
though they are not known to be finite pointwise in time for weak solutions; 
the end results may be justified in the standard way by proceeding from a 
Galerkin approximation to the solutions and then removing the regulariza- 
tion in the final results. In the usual manner 'c' and c„ are used as generic 
constants. 

Proposition 1. For Gr —t cx) and 5 > Q the Fn satisfy the differential 
inequalities 

iFo < -i^Fi + ciur^Rc^+^^Fo, (33) 

^A < ~j'^F2 + C2iy~^Fl + ci^r^Re^+^^Fi, (34) 

and, for n > 2, there are two alternative versions 

hFn < -vFn+i + c„a (ll VmIIoo + vr^Rc^+^^) F„ . (35) 

ii^„ < -ii^F„+i + c„,2 (y-^\\u\\l, + z/r2Rci+2^') Fn. (36) 

Remark 1. When the inequality F^ < F2F0 is used in H34|l . the resulting 
differential inequality for Fi demonstrates the inability of these methods, 
as they stand, to gain control over Fi for arbitrarily large initial data. 

Proof. The proof follows in four steps. 

Step 1: Let us begin with the proof of H33|l : Leray's energy inequality is 

hHo<-iyH,+Hy^\\fh (37) 

Adding and subtracting the quantity i^t^HV/IH, it is seen that 

ih < ^vF, + uT^WVnl + Hl'^Wfh (38) 

Because the forcing is narrow-band as in Tabled it is possible to reduce a 
derivative on the || V/II2 term. This, together with Young's inequality (using 
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gr" > as a parameter with g is to be suitably chosen below) to break up 
the last term 

ii^o < -^^1 + ^i?o + r2 (i.g + ^) ||/||2 (39) 

g is determined by making the coefficients of Hq and t^||/||| equal 

3 = -^ + K^"' + -"'}'^'- (40) 

With r chosen as in (j28|l with 6 > Q, g becomes 

g = r-i ([l + Gr-(2^'+i)}'^' - Gr-^'+'^A . (41) 

Consequently, g ^ t^^ as Gr ^ oo in which case 

T-^ = vI-'^gM^ < ciy£-^Re^+'^^ (42) 

In this limit (|39|l can be written as in (j^SJ- 

Step 2: The proof of H34|) is found directly from {u) = curlu) 

\Hi<-vH2+ l\j-{ui-V)udV + Hl'^\\Vf\[2. (43) 

J n 

The middle term can be estimated thus 

j uj-{uj-V)udV < ||u;|[|||Vm||2 < cF^'^fI''^ (44) 

having used the Sobolev inequality ||a;||4 < c||Va;||2 ||'*^||2 • The proce- 
dure with the forcing is then used as in Step 1 to obtain 

iFi < -uF2 + cF^'^fI'^ + cur^Rc^+^'Fr, (45) 

Young's inequality on the middle term finally gives H34() . 

Step 3: For a proof of (|35|l . consider the ladder of differential inequalities 
satisfied by the iJ„ for n > 2 (see 1^1113) 

iij„ < -uH,,+i + c„||Vm||ooH„ + Hy'||V'7||2. (46) 

(I35|l is proved by following the procedure with the forcing as in Step 1. 

Step 4: The alternative to the differential inequality (|46|l for n > 2 is 

ii7„ < -ii^iJ„+i+c„j.-i||M||^i7„ + i?,y2||v"/||2. (47) 

(I36|l is found by using the same procedure as in Step 1 except that the 
ciuantity ^i^T^|| V"+"'^/||2 is subtracted whereas vt^WV^^^ fW^ is added. 
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3.3. Long-time averages 

In f^ it was shown how the quantities k„ are ordered such that k„ < 
Kn+i- There is no k„ that is known to be a priori bounded. What is known 
is the boundedness of the long-time averages defined in |0J. The equivalent 
of Leray's bulk dissipation estimate in terms of ki instead of Fi is found 
from (|33() by dividing through Fq and long-time averaging 

e{4)<cRc'+^'. (48) 

The first of the following two theorems states results on long-time averages 
for higher values of n. This estimate can be found in |43j with a wider range 
of (5. 

Theorem 1. For Gr -^ cxd and the parameter 5 lying in the range < 6 < ^ 

3(n-l) 

71 > 1 (49) 

(50) 





£ (k„) < c„ 


{in " 


Re^ 


II 




where 


A„ is defined by 


A„ = 3 - 


5 
2^ 


S 

+ -. 
n 




Prool 


F. Step 1: Consider first K2X- 










('^2,1) = 




)<< 


/F2\ 

\F?/ 


1/2 
> 




< 


f^ V ^2 \ 


1 
21^2 


(Fi) 





(^1 



.1/2 



(51) 



where Young's inequality has been used at the last step. Dividing inequality 
(|34|l in Proposition n in H3.2l bv F^ and long-time averaging gives 

-^(||)<c.-^(FO+.r-i(Fri) (52) 

and so 

(k2,i) < c^'^ {Fi) + yr-^ (Ff 1) (53) 

The last term is 

-1 



.r-(^r^)<^5i5^^ = ^^L-3Gr3-^ (54) 



from which it is concluded that 5 must lie in the range < (5 < i to be 
certain that this term decreases as Gr — > 00. Because 

(Fi) < 1/2^3^-4^0^ (55) 

then it follows that 

^(k2,i) <c (Lri)^Re3 (56) 
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Step 2: Now consider the quantities (k„+i,„) for n > 2 



\l^n+l,n) 



Fn+1 



< 



Fn+1 



F 



2n/{2n-l) 



1/2 \ 

pl/2(2n-l) 
n 

1/2 , , 1/2 

pl/(2n-l)^ I 



(57) 



(2n-2)/(2Ti-l) j^l/(2n-l) 



> (2ri-2)/(2Ti-l) ,j-i il/(2n-l) 



Having used the fact that k„,i < k„+i,„, H57|I and (|58|l give 

-Fri+1 



(58) 



('>:n+l,n} < 



-i(2n-l)/2n 



[Z.-2 (^^)] .. (59) 



so a Holder inequality gives 

2n (.„+!.„) < (2n - 1).-^ (p^S^b) ) + '^"^ ^^^^ 



(60) 



To estimate the first long-time average on the right hand side, consider the 
second F„-ladder in l|^ 



hF^ < -^lyFn+i + c„ (iy-'\\u\\i^ + lyr 'Re) F, 



Now define 



Y _ P 2„-l 



(61) 



(62) 



and turn H61|l into a differential inequality in y„ which involves dividing 



by F, 



2n/(2n-l) 



To achieve this we use lluIlL — CK21F1 and recall that 



K2,i < Kn.i, then 



1 2TI-2 



M\lF-^"-' < CK2A K\^i] '"-' < c^ly^F,'"-' (63) 

Hence (|61fl can be rewritten as 



(n - i)(i; + lyi-^ReYn) > h.v^^^ - cv-^n^Y' f}"-^' (64) 



Making up the coefficient in v to that in H6U|I . a Holder inequality on the 
last term gives 



2 -Fr)-I-1 , s 3-2 



F, 



2n/(2n 



r„ + vr'^RcYn 



1 



+ - -{At2,i+c(2n-2)j/-2Fi}. (65) 

In — 1 
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Taking the long-time average of this in (|6()|l we have 

2n(K„+i,„) < {K2,i)+Cny~^{Fi)+y^^{2n~l)r^Y{,c{Y^) (66) 

The long-time average of y,i has vanished and the last term (F„} is bounded 
above (because F^ is bounded below) so the long-time average is zero. Thus 
when H55|) and H56|l arc used we have 

£ (k„,i) < i (k„+i,„) < c„ {Lr^f Re^ (67) 

Note that the exponents of Re and L/ (, are uniform in 7i; only the constant 
is not. (|67|l can now be used to estimate {Kn) in the final step. 

Step 3: Rewrite (k„) in the following way: 

' t)")^((t)'"<«H 



2ti-2 



2n-2 



<Ti^l)^) < (nn.i)^^ (kI)^-^ (68) 

Using our estimate for (k„.i) from (|67|l and also that for (k^) from (|48|l . 
the result in H49() is proved. 

The first infinite set of non-trivial, bounded, long-time averages were those 
found by Foias, Guillope & Temam |7]. These are related to those in Theo- 
rem^ and particularly to the estimates for k„.i in 167|l . 

Theorem 2. For Gr — > oo the long-time averaged quantities of Foias, Guil- 
lope & Temam T^ are estimated in terms of Re as 

I {\\u\\oo) < ciu {Li-^f Re"" (69) 

i/F,r^\ < cna^^^ (Lri)'Re3 (70) 

d\\yu\\]i^)<c,u'/^iLr'fRe\ (71) 



Remark 2. The case ri = 1 is distinct from the result in ^H] because of the 
(i£~^) on the right hand side. 

Proof. The proof follows from the Sobolcv inequalities 

\\u\\oo<ckIIIf^ WVuW^KcfilllFl^^ (72) 

with the estimates (|67|) for k„_i and (|32|) for (Fi). The quantities in H7U|1 
can be rewritten in terms of k„^i and Fi and the result follows. 

A lemma is now proved that will be useful in later sections: 
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Lemma 1. // any Fm (nm) is bounded on a time interval [0, T] for 1 < 
m < n then so are all Fn (nn) for n > m. 

Proof. Consider (|35|1 in Proposition Q] in H3. 21 above: for n > 3 a Sobolev 
inequality gives 

||V«|U < c ||V"«||^||V«||^-" < F^/'Ft""^^' (73) 

where a = 3/[2(n — 1)]. There is an inequaUty for the Fn of the form 

Fp"^<F^^,F^^^^ (74) 

The choice oi N = n, p — n — 1 and q = I gives 

-Fn+i<^F,r'/F{^ (75) 

so, in consequence, (p?K| becomes 

hK < ~vFt~'lFf~' + cFI+'^'^f[^'''^'^ +cvi-^Fn (76) 

Because nj(n— 1) > 1 + a/2, (|76(l makes it clear that if F\ is bounded above 
at any time then all Fn are bounded. If any Fm is bounded for m > 1 then 
Fi must also be bounded (from 1)74(1 'I . in which case all F„ are bounded 
for any n > m. The same results hold for the k„ because the divisor Fq is 
bounded from above and below. 

4. Intermittency: the binary form of the time-axis 

In the summary section, ^ it was discussed how the effective viscosity 
could be increased by proving that the ratio Kn+i/Kn has a lower bound 
that is greater than unity under certain circumstances. This was discussed 
in the context of the ladder of differential inequalities H18|l for the F„ which 
is repeated here 

/ / \2(n+l) \ 

Fn < -\y>il ^^ + CnV-^KlF^ + z^r^Rc ^„ . (77) 



2 " 



The task of this section is to investigate lower bounds on the ratio k„+i/k„. 
In the rest of this paper, the two lengths L and t will be taken such that 
i = L/2tt to reduce algebra. Additionally, the parameter (5, lying in the 
range < (5 < i, that appears in the exponents of many of the estimates 
of the previous section, will be taken as arbitrarily small (but fixed) and 
ignored hereafter. 
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Theorem 3. For the parameter fi taking any value in the range < /i < 1, 
the ratio k„+i/k„ obeys the long-time averaged inequality (n > 1) 



Hn+l 



Kn 



(LK„)^Re~-^" ) > (78) 



where the Cn are the same as those in i49l) . Hence there exists at least one 
interval of time, designated as a 'good interval', on which the inequality 



> [LKnf Rc"^" (79) 



holds. Those other parts of the time-axis on which the reverse inequality 
holds 

< [LKnf Rc"^" (80) 



are designated as 'bad intervals '. 

Remark 3. In principle, the whole time axis could be a good interval, 
whereas the positive time average in (j78|l ensures that the complete time- 
axis cannot be 'bad'. This paper is based on the worst-case supposition 
that bad intervals exist, that they could be multiple in number, and that 
the good and the bad are interspersed. This is what is meant in this paper 
by a 'potentially binary character', although the precise distribution and 
occurrence of the good/bad intervals and how they depend on n remains an 
open question. It will be left until later (Theorem EJ to prove that the bad 
intervals are finite in width. 

Proof. Take two parameters < ^ < 1 and < a < 1 such that n + a ~ 1. 
The inverses /i~^ and a^^ will be used as exponents in the Holder inequality 
on the far right hand side of 



thereby giving 

/ / Kn+1 \ \ \ ( \'^nl \ _ /, «\ / \'*n/ 



Navier-Stokes information can be injected into these formal manipulations: 
the weak solution upper bound (|49|) and the lower bound Lk„ > 1 can be 
used in the ratio on the far right hand side of H82|) to give (|78|l . with the 
same c„ as in (|49|l . 
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4-.1. Bounds within good intervals 

On the good intervals, application of the improved lower bound H79|l to 
the differential inequality H77|l appears to imply that /i must satisfy 2ii{n + 
1) > 1 for the exponent of the negative term to be larger than the positive. 
To harden this argument it is necessary to convert H77() into a differential 
inequality in Fn alone. This can be achieved because the divisor within k„, 
namely Fq, is bounded above 

+ cu-^F^F^ + i/L"2ReF„ (83) 

For n>2 arbitrarily large initial data a singularity can be prevented from 
forming if the exponent of the negative i^„-term be greater than that of the 
positive 

(1+/X)(71+1) 271 + 3 1 ,„^, 

■^ ^-^ > — ?; ^ A« > 777 7T (84) 

n 2n '^ 2{n+l) ^ ' 

as predicted. It is not possible to take the infinite time limit because of the 
finiteness of the interval but the value of Fn = Fn.max that turns the sign 
of the right hand side of H83() is bounded by 

Fn,max < L "Rc'^"^ Fq^j^gx = Ubd (85) 

_ 4n[A„(n + l) + 2] 

The exponent 7„ > when 2/i(n + 1) > 1 and we have used the fact that 
Fo^max = cLv^Ke . In terms of Figure IT^ in SjSJ it is necessary to prove 
that the solution in the good region can become large enough to form an 
initial condition for weak solutions in the bad region. This can be proved by 
the following argument: consider that on bad intervals the k„ are bounded 
below uniformly by 

[LKn{t)f > c„Re^" (87) 



where c„^^ = Cn . In terms of Fn this can be expressed as 

Fn > Cn,fj.L "Re f" Fo^min = Lbd (88) 



The question revolves around the relative sizes of the lower bound Ltd in 
(P|l and Ubd in ^ 

bd I 0,rnax \ ^-^^-^-^^^-^ ^^ Re > 1 



Lbd \ FQ,n 

Hence it is possible for Fn to reach magnitudes at the edges of the good 
region that lie above the lower bound in (|88|l . 

For the case n — 1, the following Lemma is applicable 
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Lemma 2. When n ~ \ no singularity can form on good intervals provided 

Proof. This follows immediately by applying Theorem O to H34|) . 

Nothing has yet been proved so far regarding the widths of the good and bad 
intervals, {At)b and (^i)g respectively, nor have we any further information 
regarding their nature. While it is possible that they may form pathological 
fractal subsets of the time-axis it will be assumed that these intervals are 
simple open or closed sets; the next section is devoted to estimating upper 
bounds on {At)^. Here it is shown that a lower bound can be found on the 
ratio of the average widths of the good and bad intervals. The argument 
is based on an elementary application of the Markov-Chebychev inequality. 
Consider an interval of time [ip, tq] that contains an equal number N of 
good and bad intervals of widths {At)g and {At)f^ respectively. Define the 
average widths as 

1 N -, N 



Theorem 4. Consider an interval of time [tp^ tq] containing N pairs of 
good and bad intervals. In the limits TV — > oo and [tp, tq] — > oo, provided 
{At)b > 0, the ratio (At) /{At)^ diverges as 

> c„Re^">^M-^J as Re -> cx) (91) 



m^ > c.^M^-^ 

{At) 



b 

Proof. Given H90|l . the fraction of time occupied by the bad intervals sat- 
isfies 

S-(^^)'^ <^l dt 



1 ( Lt t ]LKndt\ 

- t,-tp\ c„,^Rc^"/'^ J 

where Tp^q = [LK„(i) > c„_pRe "'^]n[<p, tq], so as A^ ^ oo and tq—tp -^ oo, 
we have 



mb . (i-«) .\,^j^A'-^ 



< — ^ — ^^-^ < 



(Z\i), + {At\ - c„,^Rc^"/'' 
where we have used (|87ll and (|49|1 . Hence we have the result. 



(93) 
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5. What happens in the bad intervals? 

It is necessary to prove that the bad intervals are of finite width: that is, 
an upper bound is required on At = t — to where io is the initial time of some 
arbitrary bad interval. Technically speaking, there should be a superscript 
label for the ith bad interval such that At = Z\i(*' and another on to = t^, 
but these have been dropped for convenience. Recall that luq = vL^^ and 

It will become necessary to solve inequalities of the type 

£{At) < uj-^Re-'^ (95) 

for /3 > as Re — > oo. It is not difficult to show that when /3 > 1, to leading 
order 

LUoiAt) < Re-^ (96) 

whereas when < /? < 1 then, to leading order 

wo(^i) < (l-/3)Re"^lnRe (97) 

The main task of this section is to show that the bad intervals have a finite 
widths and to find an upper bound on these. This requires two subsidiary 
estimates for 

e"°^'=^*Fi(i)dt and / e"«^'=^*K2,i(t) dt (98) 

At J At 

5.1. Two subsidiary estimates 

Lemma 3. An estimate for the exponentially weighted time integral of Fi 
is 

f ^uoKoAtp^^^^ ^^ < ciuLRe^ + C2v'^L-^£{At) [Rti^ + 0(Re^)] (99) 

J At 

Proof. Let us denote Ho{t) = X'^it) with Ho{to) = Xq then Leray's energy 
inequality l|ll(l for weak solutions, 

^Ho<-i^Hi+H'J^\\f\\2 (100) 

in combination with Poincarc's inequality kiHg < Hi, gives 

X<-iyklX+\\fh (101) 

Let us also denote Xf by {ki = 2'k/L) 
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which has the same dimensions as X. Integration of (|l()l(l from to to t results 
in 

X{t) < Xoe-'^'^i'^* + Xf (l - e-'^'^i'^*) (103) 

Because there is no specific knowledge of to the upper bound on H{t) is 
taken over the full time-range < t < c» which, from (|l()3|l . is 

Ho{to) < (^] Gr' < ciy^LRc^ (104) 

This is properly valid after the time when transients have died out. The 
exponential decay in H1U3|I is trivial compared to exp(a'oRc At) so we obtain 

f e'^°^'''^*X{t)dt<ciyL^/^£iAt)Re^. (105) 

J At 

Now multiply HlOOfl by e'^"^''^* and integrate by parts 

1^ [ e'^°^'"^' Hi{t) dt < [ e'^"^'"'^' (-kHo + XWfh) dt 

J At J At ^ ' 

< ^Hoito) - iifo(t)e"''^'^* + ic^oRc / i7oe""^°'^*dt 

J At 

+ ciy'^L-'^£{At)Re'^ (106) 

In the general case the negative term can be dropped, leaving 

v [ e'^°^'"^'Hi{t)dt < civ'^LKe^ 

J At 

+ C2V^L~^£{At){Rc^ + 0(Re'')) (107) 

The predominant Re -term has a correction term of 0(Re ) from the fourth 
term in H106() and another of O(Re^) from making up Hi to Fi. 

Note that the first term on the right hand side of H99|) in Lemma |3| can 
be removed if the energy has the lower bound ffo(i) ^ -ffo(^o)e^'^°^''^*: see 
fjSlfor a discussion of this. 



Lemma 4. An estimate for the exponentially weighted time integral of K2. 



1 



J At ' J At 

+ c£{At){LGr)-^Rc (108) 

Proof. The time integral of K2,i can be estimated from 

/ K2,ie^°'^'''^Ut^ f e^°'^''^\F2/F^y/^Fl/^dt 

J At ' J At 

<^ I e"°«'=^*(i^2/-Fi')dt + A / e"°^<=^*Fidt (109) 

2 J At 2zy^ J^^ 
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The first integral on the far right liand side of Hin9|) can be estimated by 
using the inequality for Fi from H34() in Proposition ^ in i )3.2l 

ih < -^F2 + ciy-^F^ + uJoRcFi (110) 

Dividing (|ll()|l by Ff , multiplying by e'^"^''^* and integrating 
-/ e''°^'"^'{F2/F^)dt<ciy-^ f e'^f^^^^Fidi 

4 J At J At 

+ i(^r'We""'"=^*-^i"'(io)) (111) 

The last term can be rewritten in terms of £{At) which leaves a F^^^ term. 
The upper bound on this is proportional to Gr^ , which can be ignored as 
small. 



J At ' J At 

as in pU8|l above 



At J At 

cL'^£{At)GT-^Rc (112) 



5.2. An estimate for {At)i, when n > 2 

The two estimates above for the weighted time integrals of Fi and ^2,1 
allow us to prove the main result of this section for n > 2. Define 



K+i {2n-2 



fj, \2n — 1 J 2n— 1 
Then a„ > 1 if /i is chosen such that 



(113) 



whereas < a„ < 1 if /^ is chosen to lie in the range^ 

f n-l\ , f 2n-2\ , , 



Theorem 5. For n > 2, if an > 1 the width of a bad interval is bounded by 

Cn^iLOoiAt)b < Re""" (116) 

whereas i/ < a,i < 1 

£n^2UJaiAt)b < Re~'^ InRe (117) 

^ Note that for n — 2 the lower bound on jj, in 111511 is greater than 21,^+1) ■ 
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Remark 4. There appears to be no obvious parallel result for the finitcness 
of bad intervals in the case n = 1, although Lcnuira [3 gives a lower bound 
/x > I for the prevention of singularities forming on good intervals. 

Proof. Let us return to Proposition ^ in H3.2I inequality (|36() . and recall 
that Wo — vL~^ 

hFn < -^vFn+i + c,, {v-^\\u\\l, + ujoRe) F„. (118) 

This was manipulated in ii3.3l to produce (|64|l which is re-stated here as 

(n - i)(r„ + L^oRer„) > jiy^^ - C21^-'k.^F^ (119) 

The first term on the right hand side of (|119|l can be estimated as 

J;S-i) ^ -^^0^ ^ cn^{.'LRer-^ (120) 

having used the fact that k„ < k„+i. This result, together with a Holder 
inequality, gives 



2n-2 
^2n-l UJqUc At 



{2^ 3 2n+l "1 

I^^^^~^K2.1 + I/"5^^^Fl M121) 



So far this has just been re-arrangement of pi8|) . The lower bound on k„ 
is now applied to the first term on the right hand side along with a time 
integration 



(n - i) {y„(i)e'^«^^^*} > c^;-l^i/5^L-i£:(Z\i)Re^57; 



-1)^ 2n~l 



2n-2 



J At 

2n-2 



J At 

-f (n-i)r„(to) (122) 

For the left hand side it is sufficient to show that this is bounded above by 
a very small number on a bad interval 

Yn = K~^F„^^ < Liy~ ^^ Re- T^^^Gt'^^ (123) 

Using Lemmas O and 01 a comparison of the major terms in H122|l shows that 






(124) 
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For n > 2 the left hand side is always positive provided /i is chosen in the 
range 

^"-^2"-2^ ' >5 (125) 



/i \2n — 1 J 2n — 1 
or in the range 

, < A„,r (1^) (126) 

To solve (|124f) use the definition of a„ in (|113|l giving 

c„wof(Z\i) < Re''^" (127) 

The solution of this depends on whether a„ lies in the range a„ > 1 or 
< a„ < 1. The estimates in (|95|l and H96|l are appropriate. 

5.3. Intersection of bad intervals: the relation to Schejfer's singular set 

Figure I2T^ of §2.4 is a representation of good and bad intervals for some 
n > 2. Since it must be assumed that the position of the intervals changes 
with n, the intersection of all the bad intervals for 7i > 2 is pertinent: only if 
this intersection is non-empty will singularities be possible. For each n > 2, 
let us designate a bad interval as the set B„ on the time-axis on which 

c„:!!!^ < iLKn)^Re-^- (128) 

Moreover, because Lk„ < Lk„_|_i, on this set there is a lower bound 

(iK„)'^>c„Re^". (129) 

Now consider the set ;B„+i on which 

c„+i^^^^ < (iK„+i)^Re'^"+i ^ (iK„+i)^ > c„+iRe^"+^ (130) 

Kn+l 

Then on the intersection Xn+i = Bn H Bn+i, we have 

ci+^c„+i(L«„+2) < (LK„)(i+^)'Re-^"+^-(i+^)^" (131) 

Using {LKn+2)'^ > (LKn+i)'^ > c„+iRe'^"+^ on X„+i, a new lower bound is 

(Lk„)1+^ > (c„+iRe^"+i) ' " (c„Re^") (132) 

Now consider the intersection X„+2 = Bn H Bn+i n Bn+2- On this set there 
is a larger lower bound 

(iK„)*'+^'' > fc„+2Re^"+0 ' " fcr,+iRe^"+0 fc„Re^") '^"^ (133) 



We wish to find a lower bound on Lk„ on the set of p intersections 

In+p = Bn n Bn+1 H . . . H Bn+p (134) 
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By inspection, the general formula for the lower bound of Lk„ on X„+p is 






c„Re^ 



c„+p+iRe^"+''+^) {nlo4l^?j Re^"'" 



where 



(135) 
(136) 






ep = Eep.: = M-M(i + A^)'+'-i} 



i=0 



K„(t) 



Re 



^n/fi 



fle^ 



1 


-4 


1 (Lk„)^ > i?e^- / 






■; / 








\ I 
\ i 

\ 1 

Long-time avaragc 




\ 


-^ 




\J 




^t 



(137) 



2=0 

We are particularly interested in the limit p ^ oo so we write 

(Lk„)'' > Re^"°°' 
where An' is defined as 

c(-'Re-"' = {(c„+P+iRe^'--+'^--'') (ilLocf.':,;)'}''^^ (138) 
and then, because Xn+p > A„ and Cn+p > Cn, it follows that 

where ^p is defined by the sum 



(139) 



(140) 



Figure 3: Similar to figure ITU representation of good/bad intervals for k„ with 
a black strip representing the bad interval used in the intersection table |^ 
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1n+S 



K.n + 7 
H„ + 6 
Kn + 5 

l^n + 3 
Hn + 2 
Kn + 1 



Table 6. The lowest continuous horizontal black strip is the bad interval of k„ 
shown as the black strip in Figure 3. The strips in next 8 levels are an illustration of 
how some randomly chosen bad intervals (k„ —r Hn+s) could appear. The thicker 
strips at the highest level are the intersection of the 9 strips below. 



The potentially singular set, 5'°°-', is given by 



5* 



Bi n ^2 n . . . n s„ n , 



(141) 



must necessarily include Bi, the singular set of ki (and therefore Fi). The 
range of values of /i expressed in (|115|) and Theorem [SJ are valid for n > 2. 



These narrow to < /i < 4 in the limit 



oo. As already pointed out in 



Lemma |5] of M.ll a corresponding separate calculation for Fi shows that /i 
lies in the range i < /i < 1 for n = 1. When the allowed ranges of /i are 
taken into account for good and bad intervals we conclude that 

Theorem 6. For all bad intervals to be finite for n > 2 and for no singu- 
larities to form in good intervals for n > 1 the allowed range of fi is 



< /i < 



(142) 



^(oo) jg related to SchefFer's potentially singular set: his set is technically the 
union of all sets S^°°^ associated with every bad interval. Scheffer showed 
that this set has zero i-dimensional Hausdorfl measure ^K], which means 
that it must consist of, at most, points. Whether the k„ actually become 
singular on this set is still an open question. From (|138|I - H14()() we have 



^Re' 



< a 



(-)Re^"°°' < 



lim 



hp+i 



Re-^ 



(143) 



Divergence in this limit would guarantee singular behaviour if the set S^°°' 
is non-empty but there is no evidence that the product in H138|l diverges in 
the limit even though the upper bound in (|143|1 is infinite. From (|143|l it is 
clear that An' > A„ so all the estimates of the previous sections dependent 
upon A„ should be replaced by An . This paper, however, furnishes no 
evidence on the distribution of the intervals; Table 6 is simply a pictorial 
representation of some randomly chosen bad intervals associated with k„ — > 
Kn+8 to illustrate how the final intersection may form. 
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5.4- Dangerous sub-intervals 

In addition to the intersection idea of the last section, we consider the 
special set of sub-intervals within each bad interval on which F„ > 0. Con- 
sider the jth sub-interval within the zth bad interval: this is designated as 
dangerous sub-interval of width {At''^-') with an initial value of designated 
as tQ-' . It is on these where singularities are possible: they are not possible 
where any one of the Fn is decreasing. (|114|l and pi5|l show that the small- 
est lower bound on A„//z is A„//i > 5. Because A„+i > A„, and replacing 
A„+i by ^„+{, define 

bn = ^^ - 4 > 1 (144) 



Theorem 7. Dangerous sub-intervals are bounded in width by (n > 2) 

uJo{Ati^) < CnRe"''" bn > 1 (145) 

and on these sub-intervals 

L-^Fi > w^Re''"+'* \\u\\oo > iu;oRe''"+'' (146) 

||Vm||oc, > c^oRe''"+'* (147) 

Remark 5. Note that 5„ > a„ so the upper bounds of these sub-intervals 
are smaller than those in Theorem 

Proof. Consider 1)16(1 

iF„ < -iiyFn+i + {c„v-'\\u\\1^+ujIRc) Fn. (148) 

Now use the Sobolev inequality ||m||^ < cHn,iFi, and divide H148|) by Fn- 
Then on these sub-intervals 

't^+i.n < CnV-'^Kn,iFi + c L-^Re (149) 
Now we know that k„+i_„ > «„+! so 

k'+i < {cnV^'^Fif + cL-2Re (150) 
Now the lower bound Lun+i > Re "+1''' is invoked giving 

'ci^^^Re^Sl/A'j' _ cRe < (cnLiy-^Fif (151) 
Because Xn/fJ- > 1 for n > 2, the Re term is small in comparison, leaving 

c(,°°'Re''"+i/^ < Liy-^Fi (152) 
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Multiplying by the exponential term, integrating over At~^ , and then using 
the exponentially time- weighted integral of Fi in Lemma 13 gives 

CnUJn£{At+) [Re^^^+i/^ - Re^j < Re'' (153) 

Now we know that yl„°^{//i > 5 and 6„ > 1 so the result in (|145|l follows. 
The definition of a„ in (|113|l guarantees that 



A 



(oo) 



bn = '-^^ - 4 > a„ (154) 

which is the correct way round. 1)146(1 and H147|) follow from (|152|1 . 



6. A conditional regularity result 

The reader who has followed the proof of Theorem |S1 will have noticed 
that attempts to prove Navier-Stokes regularity fail in the bad intervals. 
There, use was made of the variables Y„{t) = Fn ' ^ "~ •' defined in (|H^ . 
To prevent the formation of singularities, it would be necessary to show 
that Yn can never touch zero in a finite time. Within the dangerous sub- 
intervals of t i5.4l this can be achieved provided the energy is bounded below 
in a certain manner. Specifically the result is 

Theorem 8. The Navier-Stokes equations are regular if, in dangerous sub- 
intervals {At]^_^), there is a lower bound on the energy 

Ho{t) > Hoitli')e--°'"'^' (155) 

Remark 6. Over the very short time interval At the exponent on the right 
hand side of H155|l is very small, so the right hand side is almost Hq (ig ■' ) . 

Proof. They key point preventing progress regarding bounding Yn(t) away 
from zero is the set of extra terms in Lemmas|31and0]that are not coefficients 
of £{At). This creates negative terms on the right hand side of H122|) that 
cannot be controlled. To circumvent this problem it is necessary to remove 
two hurdles. The first is the last pair of terms within inequality 1)111(1 



f e"°^'=^*(i^2/i^f)di<ci/-3 /■ 

J At J At 

+ i[Fri(i)e"°^<=^*-Ffi(io)] (156) 



On dangerous sub-intervals t^-' where i*i is increasing, the last term in ((156(1 
can be re-written as 

i (Ff i(t)e-»«'^^* - Fr\tl;^)) < ic.oReFri(4^^")£(^ir) (157) 

This term is now classed as one of the £{At''_^-') terms and is merely a term 
of lower order than the dominant Re term. 
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Secondly, in (jlOBfl if it assumed that on these sub-intervals that p55|l is 
true then the extra terms in H99(l can be removed, leaving 

f g^oRe4t^^(^) ^^ < c2iy^L-^s(At\') [Rc^ + 0(Re^)] (158) 

J At 

which again is proportional to £{At''_f). Thus (|122|l becomes 

(n - 1) {y„(i)e""^'=(^*+''} >{n- hWniAf) (159) 

2,.-3 , ■ . f (2r.-2)A„ + i-4^. /| 

+ c„i/3;^L~i£:(Z\i+J) <^ Re (5^^^^^)^^ - Re^ ^ 

with no negative terms on the right hand side. Given that /i is chosen in 
the restricted ranges in (|114|l and H115() within Theorem [SJ and that £ > 
for f > 0, then Yn,{t) can never be zero. 

7. Discussion 

To summarize the arguments of this paper, it has been shown that very 
strong fluctuations in the Kn{t) can occur in time, lower bounds on which 
are much higher than the long-time average H15|l 

(Lk„) < c„Re^" A„ = 3 - — + - (160) 

2n n 

This is based on the raising of the lower bound on the ratio Kn+i/^n away 
for unity, a result which is expressed in Theorem |3 of i0 

cn^^:^ > {LK„rRe-^- (161) 

and which is effective only on the good parts of the time axis. On those 
parts of the time-axis where the reverse of (|161|) is true, no upper bound 
has been found on the k„ but very large lower bounds exist on these of the 
form 

Lk„ > c„Re^"/'' (162) 

The above results are valid for n > 2. By including intervals at n = 1 the 
intersection set of bad intervals 5'°°' is related to Scheffer's singular set of 
potential singularities, in which case the right hand side of (|162|) can be 
raised again by replacing A„ by An . The constant^ fj, is then constrained 
to the range i < /i < |. 

A picture emerges of Navier-Stokes solutions that are regular on 'most' 
of the time-axis which is punctured by short, active intervals. While no 
upper bound on k„ has yet been found within these intervals, to become 



^ In fact illblH breaks the dimensional scaling of the standard Sobolev and 
Gagliardo-Nirenberg inequalities although how the introduction of the exponent 
/i affects this in a precise manner is not yet clear. 
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singular k„ would have to find its way through a non-empty intersection 
in a similar manner as the illustration in Figure (|5.3(l . Notwithstanding the 
credence that must be given to the long-standing and widely held belief that 
the Navier-Stokes equations must be regular for arbitrarily long times, an 
alternative emerges that no upper bound exists on the k„ but the poten- 
tially singular set S^°°' in t l5.3l is either empty or only allows extremely rare 
singular events. 

The results in this paper are consistent with ideas that have existed for 
many decades concerning intermittent flows [Tlll8lf5l;,:22;,-36] but it has to be 
acknowledged that our results are lacking in four areas: 

1. The bounds are almost certainly lacking in sharpness; the two main 
places where the size of the bounds appears is in A„ and in the Re 
global estimate for Hq in H104() . which is certainly not sharp. Given 
these results, the state of the analysis is such that it may be premature 
to suggest specific numerical tests. 

2. In addition to a lack of control over the k„ within the bad intervals, 
their distribution and sensitivity to the value of n is an important but 
unanswered question. 

3. The nature of solutions in the good intervals has yet to be properly 
established. Because solutions are bounded pointwise in time and are 
also constrained by the long-time average H160|) it is to be expected that 
they should show a strong degree of quiescence, particularly within the 
central parts of these intervals. This has yet to be demonstrated. 

4. How solutions at the junctions of the good and bad intervals connect to 
each other is not clear as we have only weak solutions within the bad. 

Results of this type derived in the manner of are not confined to the 
three-dimensional Navier-Stokes equations but could be applied to simpler 
problems. All that is needed are long-time average bounds for weak solu- 
tions constructed in such a way that the equivalents of k„ are bounded 
below. An example that springs to mind is the case of the two-dimensional 
Navier-Stokes equations. Not only are these regular but tight estimates exist 
both for the attractor dimension |48j and the number of determining modes 
and nodes 021 • Other examples might be the alpha and Leray models of 
turbulence |50II51II52| or the complex Ginzburg-Landau equation |53II54| . 

The physics community has for many years used the scaling arguments 
based on Kolmogorov's original work. Frisch's book gives a detailed factual 
and historical account of these arguments |19j . Using an inertial range ar- 
gument, it has been argued in |43| that the scaling in the rigorous upper 
bounds on (k„) from ((T^ (repeated in (|16()|l above) may be interpreted in 
terms of the structure of the Fourier spectrum Eg (fc) if a scaling of the form 
Es ~ k~'^ is assumed in the inertial range, up to the cut-off wave-number 
Lkc ~ Re'^ . Disregarding the correction from S, the a priori bounds in 
(|160|l are consistent with q = 8/3 and qc = 3. Such a k~^^^ spectrum has 
arisen in at least two previous studies. Sulem & Frisch |55j have shown 
that a fc^^/'^ spectrum is the borderline steepness capable of sustaining 
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an energy cascade in the Navier-Stokes equations when the total energy is 
finite. Mandelbrot [201, and later Frisch, Suleni & Nelkin in their toy (3- 
model |56II1H| , came upon this same scaling exponent as an extreme limit of 
intermittency in the energy cascade. They found that if the energy dissipa- 
tion is assumed to be concentrated on a fractal set (in space) of dimension 
D = 8 — 3q, then the energy spectrum scaling is of the form Eg ^ fc^*. 
Within this picture, the exponent q = 8/3 thus corresponds to dissipation 
concentrated at zero-dimensional points in space. Interestingly, the conven- 
tional Kolmogorov k~^'^ spectrum for homogeneous isotropic turbulence is 
associated with D ~ 3; that is, a complete lack of intermittency with dis- 
sipation spread uniformly in space is consistent with q = 5/3. Departures 
from Kolmogorov scaling, otherwise known as anomalous scaling, can be 
associated with intermittency in the inertial-range. These arguments have 
been applied to and tested on various models such as the /3- model |56) , and 
the bifractal and multi- fractal models jT5] . While these suffer by comparison 
in not having the same degree of complexity as the Navier-Stokes equations 
- as Frisch ^Hl and Sreenivasan [^21 have both pointed out - these models 
have the merit of simplicity while capturing the main essence of the phe- 
nomena in question. More recent work on anomalous scaling has centred on 
the role of the 5*0(3) symmetry group in the expansion of the correlation 
functions |58iiM| . 
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